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A method to generate ultrahigh intense electromagnetic fields is suggested, based on the laser pulse 
compression, carrier frequency upshift and focusing by a counter-propagating breaking plasma wave, 
relativistic flying parabolic mirror. This method allows us to achieve the quantum electrodynamics 
critical field (Schwinger limit) with present day laser systems. 
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The invention of chirped pulse amplification (CPA) method and recent development of laser technology led to a 
stunning increase of the light intensity in a laser focal spot, Electrons in laser electromagnetic field become 
relativistic at intensities / ~ 10 18 W/cm 2 . The ion motion strongly affects the relativistic plasma dynamics starting 
from / > (mi/m e ) x 10 18 W/cm 2 (see Ref. |2j and references therein). Nowadays lasers produce pulses whose 
intensity is approaching to 10 22 W/cm 2 0. With further increase of intensity we shall meet novel physical processes 
such as the radiation reaction dominated regimes, which come into play at I = 10 23 - 10 24 W/cm 2 jj], and then the 
regime beyond I = 10 25 W/cm 2 where the quantum electrodynamics (QED) description is needed as the recoil of 
emitted photon momentum becomes comparable with the electron momentum 0. Near the intensity of 10 29 W/cm 2 , 
corresponding to the QED critical electric field, light can generate electron-positron pairs from vacuum 0, Q . Even 
before that limit, vacuum begins to act nonlinear ly such as vacuum polarization. These nonlinear effect have attracted 
a great deal of attention since Q: they lie outside the scope of perturbation theory and shed light on the nonlinear 
quantum electrodynamics properties of the vacuum. There are several ways to achieve such an intensity. One way was 
demonstrated in the experiments , where high-energy bunch of electrons interacts with counterpropagating intense 
laser pulse. In the reference frame of electrons the electric field magnitude of the incident radiation was approximately 
25% of the QED critical field. 

A technically feasible way is to increase the power of the contemporary laser system by some 7 orders of magnitude 
through megajoule lasers albeit quite expensive. Another way is to increase the frequency of the laser radiation and 
then focus it onto a tiny region. In this method X-ray lasers can be used To achieve more "moderate" intensities, 
10 24 — 10 25 W/cm 2 , another scheme was suggested in Ref. where a quasi-soliton wave between two foils is pumped 
by the external laser field up to an ultrahigh magnitude. Another method is based on the simultaneous laser frequency 
upshifting and the pulse compression. These two phenomena were demonstrated in a broad variety of configurations, 
where they were caused, in general, by different mechanisms. In particular, the wave amplification reflected at the 
moving relativistic electron slab was discussed in Refs. jjjj (based on the frequency up-shift of radiation reflected at 
the relativistic mirror, as predicted by A. Einstein in Ref. H3)j the backward Thompson scattering at relativistic 
electron bunch was considered in Refs. |l3(: the reflection at the moving ionization fronts has been studied in Refs. 

14]; "photon acceleration" schemes with co-propagating laser pulses in underdense plasma were examined in Refs. 

15j : various schemes of the counter-propagating laser pulses and the use of parametric amplification process were 
discussed in Refs. [Tfi| . 

In the present paper we consider a plasma wakefield in the wave-breaking regime as a tool for generating a coherent 
radiation of ultra-high intensity. Compared to the previously discussed schemes this regime demonstrates both the 
robustness and coherence of the transformed laser light. 

We examine the following scenario. A short intense laser pulse (the 11 driver pulse") induces wakefield in a plasma. 
As it is well known [17J, the wakefield phase velocity v p h = P p hC equals the laser pulse group velocity, which is close to 
the speed of light in a vacuum when the laser pulse propagates in the underdense plasma. The corresponding Lorentz 
factor is j p h = (1 — Z^) -1 ' 2 ~ ^>d/^>pe, where u>d is the driver pulse frequency, ui pe is the Langmuir frequency. The 
nonlinearity of strong wakefield causes a nonlinear wave profile, including a steepening of the wave and formation 
of localized maximums in the electron density - the spikes [18). This amounts to wavebreaking regime (see Ref. 
2| and references therein). Theoretically the electron density in the spike tends to infinity, but remains integrable 
2j. Sufficiently weak counter-propagating laser pulse (the "source pulse") will be partially reflected from the density 
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maximum. The reflection coefficient scales as "f p h and the reflected wave vector-potential scales as 7~/f, as it is shown 
below. As we see, the electron density maximum acts as a mirror flying with the relativistic velocity v p h w c. The 
frequency of the reflected radiation is up-shifted by factor (1 + f3 pll )/ (1 — (3 ph ) « &Y ph , in accordance with the Einstein 
formula [12j . It is important that the relativistic dependence of the Langmuir frequency on the driver pulse amplitude 
causes parabolic bending of constant phase surface of the plasma wave, since the driver pulse has a finite transverse 
size [l9|. As a result, the surface where the electron density is maximal has a shape close to a paraboloid. Because 
we have a curved mirror, the frequency uj s of the reflected radiation depends on the angle: 

;Us, (1) 
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where uj s is the source pulse frequency and 9 is the angle between the reflected wave vector and the direction of the 
driver pulse propagation in the laboratory frame. The curved mirror focuses the reflected light. The focal spot size is 
of the order of the diffraction limited size. In the reference frame of the wakefield it is A' s = A s ((l — P p h) / (1+ Pph)) 1 ^ 2 ~ 
X s /2"f p h, where A s is the wavelength of the source pulse. In the laboratory frame the focal spot size is approximately 
X s /Ay^ h along the paraboloid axis, and rj X s /2^ p f l in the transverse direction. In the focal spot the resulting intensity 
gain factor scales with j p h as 7^ x {lj s /uj s ) 2 x (D s /X' s ) 2 = 64(D S / X s ) 2 jp h , where D s is the diameter of the efficiently 
reflected portion of the source pulse beam. This value can be great enough to substantially increase the intensity of 
the reflected light in the focus, even up to the QED critical electric field. 

In order to calculate the reflection coefficient, we consider the interaction of an electromagnetic wave with a 
maximum of the electron density formed in a breaking Langmuir wave. In the laboratory frame, this interaction can 
be described by the wave equation 



O tt A z -c*AA z + 4 ™ n ( X - V ^ A z = 0, (2) 
m e j e 

where A z is the z-component of the vector potential, j e is the electron Lorentz factor, and 7 e w j p h near the maximum 
of the density in the wakewave wavebreaking regime. 

According to the continuity equation dtn e + div(n e t> e ) = 0, the electron density in the stationary Langmuir wave 
is given by n — n$v p h/ {v p h — v e ), where the electron velocity v e varies from — v p h to v p h (see Ref. 18j), and the 
electron density varies from the minimal value = no/2 to infinity (intcgrable) . For the breaking plasma wakewave, 
in every wave period approximately a half of electrons are located in the spike of the electron density. Therefore we 
can approximate the electron density by n (x — v p ht) — (1 + X p S(x — v p i l t))no/2, where X p is the wakefield wavelength 
and S(x) is the Dirac delta function. This approximation is valid when the density maximum thickness is sufficiently 
less than the collisionless skin depth c/u> pe and source pulse wavelength in the wakefield rest frame, i. e. when the 
wakefield is close to the wave-breaking regime. 

In the reference frame comoving with the plasma wakewave, Eq. (0) has the same form. The Lorentz transformation 
to this frame is given by t' = (t — v p hx/c 2 )j p h, x' = (x — v p ht)~f p h, y' = y, z' = z. 

We seek for solution to Eq. J5J in the form A z — A(x') exp {i(ui' s t' + k' x x' + k' y y' + k' z z')) , where ui' s — (lu s + 
v p hk x )~1ph, k' x = (k x + v p hLo/c 2 )j p h, k'j_ = kj_ are the frequency and wavevector in the moving frame, and k' x > 0. 
Using this ansatz, from Eq. Q in the moving frame we obtain 



d 2 A 



dx 



12 



q 2 A = X S(x')A 7 (3) 



where q 2 — us 12 /c 2 — k' ± 2 — uj 2 e / (2c 2 ^ p h) > and \ — w p e A P /c 2 . This equation is equivalent to the scattering problem 
at the delta potential. The solution is A(x') = exp(iqx') + p(q) exp(-iqx') for x' > (incident and reflected wave), 
and A{x') — r(q) exp(iqx') for x' < (transmitted wave), where p(q) = — W (\ + 2ig) and r(q) — iqj (x + 2iq). In a 
nonlinear Langmuir wave, its wavelength depends on the wave amplitude |18| . and for the breaking wakewave we have 
X p w ^{2 r ) p } l ) 1 / 2 c/ui pe . In this case \ = 4(27 p / I ) 1 / 2 w pe /c. Taking lj' 2 — ij^w 2 into account, we find that the reflection 

coefficient, defined as a ratio of the reflected to the incident energy flux, in the co-moving frame is « (ujd/io s ) 2 /2^ h . 
In the laboratory frame it is 

R w 8 lph {uja/ujs) 2 . (4) 

The intensity I s f in the focal spot of the source pulse, reflected and focused by the electron density maximum in 
the laboratory frame, is increased by the factor of the order of 

I sf /I s w 32(u Jd /u Js ) 2 (D s /X s ) 2 j 3 ph . (5) 
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FIG. 1: The electron density in the wake of the driver laser pulse at t = 14 x 2ir/uid- The (x, y = — 6A, z)-plane: density profile 
along the symmetry axis. Blue curves for density values n = 0.12, 0.24, 0.36n cr on the corresponding perpendicular planes of 
symmetry. Isosurfaces for value n = 0.15n cr , "blue gas" for lower values. 

Theoretically, the actual gain can be even greater, because a) the estimation (JIJ corresponds to one-dimensinal case, 
whereas the density modulation in the 3D breaking wakewave is stronger, b) the reflectance Q of the 3D paraboloidal 
mirror is greater at the periphery. 

We consider the following example. A one- micron laser pulse {driver) generates wakefield in a plasma with density 
n e = 10 17 cm -3 . The corresponding plasma wavelength is X p w 100/im. The Lorentz factor, associated with the phase 
velocity of the wakefield, is estimated as j p h ~ LUd/uj pe ~ 100. The counter-propagating one-micron laser pulse with 
intensity I s = 10 17 W/cm 2 (source) is partially reflected and focused by the wakefield cusp. If the efficiently reflected 
beam diameter is D s = 200/im, then, according to Eq. the final intensity in the focal spot is I s f w 1.5xl0 29 W/cm 2 . 
The driver pulse intensity should be sufficiently high and its beam diameter should be enough to give such a wide 
mirror, assume Id = 10 18 W/cm 2 and Dd — 800//m. Thus, if both the driver and source are one- wavelength pulses, 
they carry 17J and 0.1J, respectively. We see that in an optimistic scenario the QED critical electric field may be 
achieved with the present-day laser technology! 

To demonstrate the feasibility of the effect of the light reflection and focusing by the breaking wakewave, we 
performed three-dimensional particle-in-cell (PIC) simulations using the code REMP (Relativistic Electro-Magnetic 
Particle-mesh code) based on scheme 20]. In the simulations the driver pulse propagates in the direction of the x-axis. 
Its dimensionless amplitude is ad — 1.7 which corresponds to peak intensity 4 x 10 18 W/cm 2 x (l/im/A^) 2 , where Ad is 
the driver wavelength. The driver is linearly polarized along the z-axis, it has the gaussian shape, its FWHM size is 
3Ad x 6Ad x QXd- The source pulse propagates in the opposite direction. Its wavelength is two times greater than the 
driver wavelength, A s = 2Ad- The source pulse amplitude is chosen to be small, a s = 0.05, to reduce the distortion of 
the wakewave. The pulse shape is rectangular in the x-direction and Gaussian in the transverse direction, its size is 
6Xd x 6 Ad x 6 Ad. To distinguish the electromagnetic radiation of the driver from the source pulses, we set the source 
pulse to be linearly polarized in the direction perpendicular to the driver polarization, i. e. along the y-axis. The 
laser pulses propagate in the underdense plasma slab with the electron density n e = 0.09n cr , which corresponds to 
the Langmuir frequency w pe = 0.3ud- The plasma slab is localized at 2Xd < x < 13Ad in the simulation box with size 
22Xd x 19.5Ad x 19.2Ad- The simulations were carried out on 720 processors of the supercomputer HP Alpha Server 
SC ES40 at JAERI Kansai. The mesh size is dx = Ad/100, total number of quasiparticles is 10 10 (ten billion). The 
boundary conditions are absorbing on the x-axis and periodic in the transverse directon, both for the electromagnetic 
fields and quasi-particles. We emphasize that the simulation grid must be and in fact was chosen to be fine enough 
to resolve the huge frequency up-shift given by Eq.QJ, exhausting all the supercomputer resources. 

The simulation results are presented in Figs. ^a-ndE! Fig. ^ shows the plasma wakewave induced by the driver 
laser pulse as modulations in the electron density. We see the electron density cusps in the form of paraboloids. They 
move with velocity v p h ~ 0.87c, the corresponding gamma-factor is j p h ~ 2. Their transverse size is much larger 
than the wavelength of the counterpropagating source pulse in the reference frame of the wakefield. As seen from the 
electron density profile along the axis of the driver pulse propagation, the wakewave dynamics is close to wave-beaking 
regime. Each electron density maximum forms a semi-transparent parabolic mirror, which reflects a part of the source 
pulse radiation. 

In Fig. Elwe present the electric field components. The driver pulse is seen in the cross-section of the z-component 




FIG. 2: The cross-sections of the electric field components. The (x,y,z — 0)-plane: E z (x,y,z — 0) (green-brown colorscale), 
the plane (x, y = 0, z): E y (x, y = 0, z) (blue-red colorscale) at t — 16, 18, 20, 22 x 2n/ujd (top-down). 

of the electric field in the (x, y,z — 0)-plane. The source pulse and its reflection is seen in the cross-section of the 
y-component of the electric field in the (x, y — 0, z)-plane. The part of the source pulse radiation is reflected from 
the flying paraboloidal mirrors, then it focuses yielding the peak intensity in the focal spot, and finally it defocuses 
and propagates as a spherical short wave train, whose frequency depend on the wave vector direction, in agreement 
with Eq.QJ. This process is clearly seen in the animations produced from the data (see authors' website). The main 
part of the reflected light power is concentrated whithin the angle ~ hence this coherent high-frequency beam 

resembles a searchlight. The reflected part has the same number of cycles as the source pulse, as expected, since 
it is Lorentz invriant. The wavelength and duration of the reflected pulse are approximately 14 times less than the 
wavelength and duration of the source pulse, in agreement with Eq.Q since (1 + (3 p h)/(l — (3 p h) ~ 14.4. The focal 
spot size of the reflected radiation is much smaller than the wavelength of the source pulse. The electric field in the 
focal spot is approximately 16 times higher than in the source pulse. Therefore, the intensity increases 256 times in 
agreement with estimation (J5J. 

We emphasize that the efficient reflection is achievable only when the wakcfield is close to the wave-breaking regime 
and the cusps in the electron density are formed. As we see in the simulations, the reflection and focusing is robust 
and even distorted (to some extent) wakewave can efficiently reflect and focus the source pulse radiation. We also 
observe that despite the moderate reflection coefficient, the colossal frequency up-shift and focusing by a sufficiently 
wide (transversely) wakewave give us a huge increase of the light intensity. 

Similar processes may occur in laser-plasma interation spontaneously, e.g. when a short laser pulse exciting plasma 
wakewave is a subject of the stimulated backward Raman scattering or a portion of the pulse is reflected back from 
the plasma inhomogeneity. Then the backward scattered electromagnetic wave interacts with plasma density modu- 
lations in the wakewave moving with relativistic velocity. According to scenario described above, the electromagnetic 
radiation, reflected by the wakewave, propagates in the forward direction as a high-frequency strongly collimatcd 
(within the angle ~ l/7 P h) electromagnetic beam. 

We have proposed the scheme of the relativistic plasma wake caustic light intensification, which can be achieved 
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due to the reflection and focusing of light from the maximum of the electron density in the plasma wakewave at 
close to the wave-breaking regime. The presented results of 3D PIC simulations provide us a proof of principle of the 
electromagnetic field intensification during reflection of the laser radiation at the flying paraboloidal relativistic mirrors 
in the plasma wakewave. With the ideal realization of the described scheme we can achieve extremely high electric 
fields (in the laboratory reference frame) aproaching the QED critical field with the present-day laser technology. We 
envision the present example is just one manifestation of what we foresee as the emergence of relativistic engineering. 

We appreciate the help of APRC computer group. We thank V. N. Bayer, M. Borghcsi, J. Koga, K. Mima, G. 
Mourou, N. B. Narozhny, K. Nishihara, V. S. Popov, A. Ringwald, V. I. Ritus, V. I. Telnov, and M. Yamagiwa for 
discussion. 



* Also at General Physics Institute RAS, Vavilov str. 38, Moscow 119991, Russia 

' Also at Moscow Institute of Physics and Technology, Institutskij per. 9, Dolgoprudny Moscow region 141700, Russia; 
Electronic address: timur@apr.jaeri.go.jp 
[1] G. A. Mourou, C. P. J. Barty, and M. D. Perry, Phys. Today 51, No. 1, 22 (1998). 

[2] S. V. Bulanov, et. al. in: Reviews of Plasma Physics. Vol. 22, ed. by V.D. Shafranov (Kluwer Academic / Plenum Publishers, 

NY, 2001), p. 227. 
[3] A. Zhidkov, et al., Phys. Rev. Lett. 88, 185002 (2002). 
[4] S. V. Bulanov, T. Zh. Esirkepov, J. Koga, and T. Tajima, in preparation. 

[5] W. Heisenberg and H. Z. Euler, Z. Phys. 98, 714 (1936); J. Schwinger, Phys. Rev. 82, 664 (1951). 

[6] E. Brezin and C. Itzykson, Phys. Rev. D 2, 1191 (1970); N. B. Narozhny and A. I. Nikishov, Sov. Phys. JETP 38, 427 ( 

1974); V. S. Popov, JETP 94, 1057 (2002). 
[7] C. Bula, et al., Phys. Rev. Lett. 76, 3116-3119 (1996). 
[8] T. Tajima and G. Mourou, Phys. Rev. ST Accel. Beams 5, 031301 (2002). 

[9] A. Ringwald, Phys. Lett. B 510, 107 (2001); R. Alkofer, et al., Phys. Rev Lett. 87, 193902 (2001); T. Tajima, Plasma 
Phys. Rep. 29, 207 (2003). 
[10] B. Shen and M. Y. Yu, Phys. Rev. Lett. 89, 275004 (2002). 

[11] K. Landecker, Phys. Rev. 86, 852 (1952); L. A. Ostrovskii, Soviet Physics Uspekhi 116, 315 (1976). 
[12] A. Einstein, Ann. Physik 17, 891 (1905). 

[13] F. R. Arutyunian and V. A. Tumanian, Phys. Lett. 4, 176 (1963); Y. Li, et. al., Phys. Rev. ST Accel. Beams 5, 044701 
(2002). 

[14] V. I. Semenova, Sov. Radiophys. Quantum Electron. 10, 599 (1967); W. B. Mori, Phys. Rev. A 44, 5118 (1991); R. L. 

Savage, Jr., et al., Phys. Rev. Lett. 68, 946 (1992). 
[15] S. C. Wilks, et al., Phys. Rev. Lett. 62, 2600 (1989); C. W. Siders, et al., Phys. Rev. Lett. 76, 3570 (1996); Z.-M. Sheng, 

et al., Phys. Rev. E 62, 7258 (2000). 
[16] G. Shvets, et. al., Phys. Rev. Lett. 81, 4879 (1998); Y. Ping, et. al., Phys. Rev. E 62, R4532 (2000); P. Zhang, et al., Phys. 

Plasmas 10, 2093 (2003); N. J. Fisch and V. M. Malkin, Phys. Plasmas 10, 2056 (2003). 
[17] T. Tajima and J. Dawson, Phys. Rev. Lett. 43, 262 (1979). 
[18] A. I. Akhiezer and R. V. Polovin, Sov Phys. JETP 30, 915 (1956). 
[19] S. V. Bulanov and A. S. Sakharov, JETP Lett., 54, 203 (1991). 
[20] T. Zh. Esirkepov, Comput. Phys. Comm. 135, 144 (2001). 



